Non-Universal Quantities from Dual RG Transformations 
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Using a simplified version of the RG transformation of Dyson's hierarchical model, we show that 
one can calculate all the non-universal quantities entering into the scaling laws by combining an 
expansion about the high-temperature fixed point with a dual expansion about the critical point. 
The magnetic susceptibility is expressed in terms of two dual quantities transforming covariantly 
under an RG transformation and has a smooth behavior in the high-temperature limit. Using the 
analogy with hamiltonian mechanics, the simplified example discussed here is similar to the anhar- 
monic oscillator, while more realistic examples can be thought of as coupled oscillators, allowing 
resonance phenomena. 
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One important contribution of the Renormalization 
Group (RG) method is to show that there exists a close 
connection [y between statistical mechanics near critical- 
ity and Euchdean field theory in the large cut-off limit. 
In this approach, the determination of the renormalized 
quantities at zero momentum amounts to the determina- 
tion of a certain number of parameters appearing in the 
scaling laws. Some of these parameters are universal (the 
critical exponents) and much effort has been successfully 
devoted to their calculation. On the other hand, new 
techniques need to be developed in order to reliably cal- 
culate the non-universal parameters. 

We limit here the discussion to the case of scalar field 
theories with a lattice regularization (spin models). This 
class of models has several important applications in 
particle physics (e.g., self-interactions in the Higgs sec- 
tor) and condensed matter physics (e.g., ferromagnetism) 
which require an accurate non-perturbative treatment. 
For /3, the inverse temperature (or the hopping param- 
eter), close to its critical value /3c, one can express the 
magnetic susceptibility (zero-momentum two-point func- 
tion) with an expression which, when D < A, takes the 
form (I 



X~ (/3c-/3)"^(Ao+^i(/3c-/?)^ + ...) 



(1) 



where the non-universal quantities , , . . . are func- 
tions of the other ( "bare" ) parameters of the theory. Fol- 
lowing the discussions of Ref. , we can use Eq. (0) to 
obtain a non-perturbative definition of the renormalized 
mass Tod of the form 



A^ + A,{^y 



(2) 



for a scale of reference Kr, and a UV cut-off A. Similar 
considerations apply to the other renormalized quanti- 
ties which can be obtained from the higher point func- 
tions. In order to complete in a quantitative way this 
non-perturbative renormalization program, one needs to 
be able to calculate the non-universal quantities in Eq. 
(0) as well as the universal ones. 



This task can be achieved [01 in the case of the well- 
studied hierarchical model 10-01 ■ Using the numerical 
methods developed in Ref. one can calculate the sus- 
ceptibility at various values of (3 and extract the unknown 
parameters in Eq. (|^) by direct fits j|] . This is a rather 
tedious procedure involving successive numerical refine- 
ments. A more satisfactory approach consists in expand- 
ing about the fixed point calculated by Koch and Wittwer 
Q . In a system of coordinates where the fixed point is 
at the origin and the axes coincide with the eigenvectors 
of the linearized transformation, the RG transformation 
reads 
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where the Am are the eigenvalues of the linearized RG 
transformation (which yield the critical exponents) and 
the are calculable coefficients. In Ref. we found 
that the direct fits and the linearization agree with 12 
significant digits for the leading exponent 7. The lin- 
earization method does not provide a way to calculate 
the non- universal quantities (Aq, Ai . . .). It would be a 
great accomplishment to show that these non-universal 
quantities could be calculated accurately by taking into 
account the non-linear terms in Eq. (||). In this Letter 
we show that such a calculation can be performed in a 
one-variable version of Eq. (H) which is justified in the 
next paragraph. Furthermore, some of the calculations 
can be performed much more efficiently by combining the 
expansion described above with a dual expansion which 
can be identified with the high-temperature expansion. 

In the following, we consider a recursion relation for 
the magnetic susceptibility which reads 
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in order to approximate 



(4) 



where c = 

Z?— dimensional model and n stands for the fact that the 
susceptibility is calculated with a number of sites 2". In 
the following, we limit ourselves to a range of parame- 
ters corresponding to ferromagnetic interactions in the 
symmetric phase, and such that an infinite volume limit 



1 



exists. This means < /? < /?c (the value Pc is calcu- 
lated below) and < c < 2. Eq. can be obtained 
as follows. First, we consider the recursion formula for 
the hierarchical model in the approximation where the 
Fourier transform of the local measure is approximated 
by a polynomial of degree 2 (this is called Imax = 1 in Rcf. 
0) and then we expand the resulting recursion formula 
for the susceptibility to first order in (3. The variable is 
then rescaled in order to obtain a recursion formula in 
terms of the physical quantity x- The recursion formula 
Eq. (^ becomes an accurate approximation of the ex- 
act recursion formula for Dyson's model when n is large 
enough. A related formula is used in Ref. to estimate 
the finite volume effects (see Eq. (5.1) therein). In the 
following, we use the notation ^ for c/2. For definiteness, 
we will take the initial value xo = 1- 

The explicit dependence on n and (3 in Eq. (^ can 
be eliminated by introducing /i„ = a^"x„. The constant 
of proportionality a can be fixed by requiring that the 
fixed points of the the KG transformation in terms of the 
new variable are and 1. This yields a = /3c^/(8(2 — c)). 
The initial value Hq = 1 (the unstable fixed point), corre- 
sponds to the choice j3 — 8(2 — c)/c^. In summary 



(5) 



and the recursion formula then becomes a simple 
quadratic map (called the "/i-map" hereafter) 



hn+l = iha + (1 - i)hl , 



(6) 



together with the initial condition = /3//?c- The 
restriction to < /3 < /3c corresponds to the range 
< /iQ < 1 which implies that for positive and finite 
n, hn stays within this interval. Note that Eq. (||) 
can be used to give /i„ as a function of hn+i- This 
quadratic equation has two solutions; however, if we re- 
quire < < 1, only one solution is acceptable 
and a unique inverse can be obtained by this restriction. 
If we impose this restriction, the term "group" in RG can 
be understood in its proper sense. 

We now discuss the two fixed points. The fixed point 
/iQ = corresponds to the choice (3 = and is called the 
high-temperature (HT) fixed point. Remembering that 
< ^ < 1, we see that the HT fixed point is stable, 
with eigenvalue ^ in the linear approximation. Using a 
graphical representation of the quadratic map, one sees 
that the HT fixed point is globally attractive for the in- 
terval (0,1). At the other end, the fixed point = 1 
corresponds to the choice /3 = (3c, and is called the crit- 
ical point. This fixed point is unstable, with eigenvalue 
A = (2 — ^). Note that if ^ is fixed by our initial choice 
of the dimensionality parameter D, the value of A can 
be seen as an approximate value for largest eigenvalue Ai 
of the hierarchical model. This value is not too far off 
numerically. For instance for D = 3,2 — ^ — 1.37 which 
can be compared with the known |^J^ value 1.42717 .... 

We can expand the h-map about the unstable fixed 
point by using the reparametrization /i„ = 1 — (i„. Note 



that do = {(3c — (3) / (3c is the variable which appears in the 
parametrization of the susceptibility given by Eq. (|^). 
The recursion formula for d„ reads: 



dn+l = + (1 - X)dl , 



(7) 



with A = 2 — ^. We will call this map the "d-map". 
This map can be seen as a one variable version of Eq. 
(^. Note the similarity with the original h-map. One 
can introduce a duality relation between the two maps 
which interchanges /i„ ^ d„ and ^ 4-> A. In the follow- 
ing, we use the notations /i„ = dn and ^ = A to express 
the quantities appearing in the h-map as dual to the one 
appearing in the d-map. If the duality transformation is 
applied twice, one returns to the original quantities. 

Recalling that < ft-o < 1, we also have < do < 1 
with small values (approaching from above) in one vari- 
able corresponding to "large" values (approaching 1 from 
below) values in the dual variable. We would like to con- 
struct an expression for x which is accurate for both small 
and large values of do. In order to do this, we need to 
use Eq. (^ beyond the linear approximation. In the 
linear approximation, which is justified when do is very 
small (/? close to Pc), dn — A"do. The Hnear approxi- 
mation breaks down for values of n of order n* defined 
by the relation A" do = 1. For n larger than n*, the 
non-linear terms become important and d„ approaches 
1 from below as dictated by the global attractiveness of 
the HT fixed point. For n large enough, the linearized 
/i— map can be used to show that the HT fixed point is 
reached exponentially fast. The linearization about the 
critical point provides the usual type of expression for 
the critical exponent 7: since x j 



7--lne/lnA. 



(8) 



In order to refine the order of magnitude estimate given 
by the leading singularity, we will express d„ as a func- 
tion of do. For this purpose, we first construct a function 
y{d) which transforms covariantly under Eq. (R) : 



yiXd+il-X)d') = Xyid) 



(9) 



If we add the requirement that for small values d, y{d) ~ 
d, this equation has a unique solution as a power series 
in d: 
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A(A + 1) A2(l + A)(l + A + A2) 



(10) 



The inverse function can be constructed similarly. In 
both cases, the coefficients can be calculated by simple 
recursion relations, easily implementable on a computer. 
We can now write: 



rfn=y"'(A"2/(do)) 



(11) 



The idea of using intermediate variables with simple 
transformation properties has a long history, for instance 
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the angle-action variables in hamiltonian mechanics and 
the normal form of differential equations appearing in 
Poincare's dissertation. For continuous RG transforma- 
tions, Wegner |^ introduced the notion of scaling vari- 
ables. In the case discussed here, the construction of y 
can be seen as a discrete version of Wegner's procedure. 

Much can be said about the convergence properties of 
y{d) and its inverse. A numerical analysis of the coeffi- 
cients indicates very clearly that is an entire func- 
tion, while y is analytical on the open disk of radius 1 
and has a power singularity when its argument tends to 
1. Consequently, when < c?o < Ij one can always find 
an accurate expression for dn by using sufficiently many 
terms in the expansions of y and y^^. Note that Eq. 
( |ll| ) can also be used at negative values of n, providing 
the inverse transformations, which can be uniquely de- 
fined by requiring - as in the case of the /i-map discussed 
above - that the preimage lies in the (0,1) interval. With 
this requirement, the transformation of Eq. (^ becomes 
a group and the function y{d) a non-unitary represen- 
tation of this group. Note that one could also define a 
continuous transformation by extending n to all the real 
values. 

Everything we have done for the d-map can be repeated 
almost verbatim for the /i-map. We can construct a dual 
function y{d) transforming covariantly under the /i-map 
with an expansion in d of the form of the one given in Eq. 
( p^ ) but with A replaced by A = ^. As for y{d), we have 
clear numerical indication that y{d) is analytical on the 
open disk of radius 1 and has a power singularity when 
its argument tends to 1. On the other hand, its inverse 
is not entire but has a finite radius of convergence with 
a square root behavior at the intersection of the bound- 
ary of the disk of convergence and the negative real axis. 
Recalling that < ^ < 1 and that do = = /^//3c we see 
that for n large enough, we can use the above-described 
expansions to calculate 



In the limit where n becomes infinite, the argument of 
y~^ goes to zero and it is justifiable to retain only the first 
term of its expansion. Using the definition of the suscep- 
tibility of Eq. (H), we find that the ^-dependence cancels 
and that we obtain the high-temperature expansion: 



X 
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= 1 + /3 



4(2 - c) 



(13) 



This expansion has features which are in qualitative 
agreement with the actual HT series of the hier- 
archical model. 

One can in principle use this HT expansion to extract 
the leading and subleading singularities of x- However, 
this procedure is in general very inefficient because small 
physical effects can be amplified dramatically in this ex- 
pansion. For instance for A ~ 1.8, the sequence of ratios 
of successive coefficients is completely "noisy" and no in- 
formation can be extracted from it. When A is lowered. 



the "noise" decreases and takes the form of smooth log- 
periodic oscillating terms as in the examples discussed in 
Ref. §. 

Instead of using the HT expansion, we would like to 
have an expansion in terms the dual variable do- Such 
a goal can be achieved by combining the two covariant 
quantities y and y into one invariant quantity which we 
call A below. Using the definition of 7 given in Eq. (||) , 
one sees that A"*" = A^^ and consequently 



A = (yidnWyiK) 



(14) 



is n- independent. A can be called a constant of motion 
or an RG invariant. We can now rewrite 



A 



X 



(l-do)(y(do))^ 



(15) 



We will show later that A is a bounded function for 
< do < 1- Eq. (|l^) makes us suspect that x has 
a singularity when do becomes close to 1, or in other 
words, when (3 becomes small. On the other hand, we 
know that in this limit x = 1. This apparent difhculty 
can be resolved by noticing that y has a singularity with 
power —7 = In A/In A, consequently the dual quantity y 
has a power singularity with dual exponent: 7 = I/7. 
Consequently, {y{do))'' at the denominator cancels the 
singularity and the expansion extends globally. 

We now calculate A expressed as a function of ?;(do) = 
yQ. The invariance of A under a RG transformation im- 
plies the discrete scale invariance: 



A{Xyo) = A{yo) , 
and the Fourier mode expansion: 

4- 00 



Vo 



(12) 

^ ' with uj — 27r/lnA, and consequently 



I r^Va 



'A{yo) 



(16) 



(17) 



(18) 



The lower value ya of the integration interval is arbitrary 
and we we can choose it at our convenience and construct 
a decent series approximation for A which is accurate in 
the integration interval. More explicitly, we can rewrite 



Aiyo) = iyorm~y-\yo))) 



(19) 



and use the series expansions for y~^ and y. If yo is small, 
we need a few terms for y~^ and many for y. If yo is large, 
we need many terms for y~^ and a few terms for y. These 
two extreme possibilities are very inefficient ways to cal- 
culate the Fourier coefficients. We have compared the 
approximate values ao(rn,m) obtained from expansions 
of Eq. ( p^ with m terms for y~^ and m for y with an 
accurate value of ao and found the approximate behavior 



3 



|ao(rTT., m) — OqI oc exj>{—Ki{m + rh) + K2{m — m) ) , 



(20) 



where Ki and K2 are positive constants. Eq. ( po| ) im- 
plies that for m + m fixed, it is very advantageous to pick 
the "self-dual" option m ~ to. 

For values of ^ not too small, the contribution of the 
non-zero Fourier modes to the susceptibility is exponen- 
tially suppressed. We found indications for the following 
behavior: 



a„|/|ao| oc exp(-|n|7rw/2) , 



(21) 



as in another example of function with log-periodic os- 
cillations discussed in section 5 of Ref . |^ . The first in- 
dication is the shape of the basin of attraction of the 
stable fixed point of the /i-map, in the complex /i-plane. 
Near the unstable fixed point (/i = 1, d = 0), we can 
linearize y ~ d = 1 — /i. From Figure |l|, we see that if 
d — —5exp{id) with < S << 1, the points such that 
6*1 < tt/2 are attracted to zero. Given the behavior of 
yinuj jj-^ ^Yiis region, this requires for large values of n the 
suppression given by Eq. (21). We have also checked ex- 
plicitly for n = I and 2, that this exponential suppression 
provides a good fit of the data for w > 11. Despite the 
suppression, the non-zero modes are quite visible in the 
HT expansion because of a factor (r(7 -I- inuj))~^ which 
appears in the the expression of the HT coefficients (see 
Eq. (3.7) of Ref. [^) and cancels, in leading order, the 
suppression from Eq. (pi|). 




Re 

FIG. 1. Boundary of the basin of attraction of for the 
complex /i-map with ^=0.5. 

The construction of covariant quantities extends easily 
to the case of several variables and can be used to at- 
tempt accurate calculations with the hierarchical model. 
The extension to the case of nearest neighbor model is a 
more difficult procedure. It is clear that one can as in Ref. 
[0 use polynomial approximations for the Fourier trans- 
form of the local measure and write a recursion formula 
for a block-spin transformation. However this procedure 
generates non-local interactions. It is not clear if these 
interactions can be exponentiated in a compact way (as 
in the Gaussian case). The present work should be seen 
as an encouragement to attack this difficult question. 



Two remarks can be made regarding the construction 
of covariant quantities in the multivariable case. The 
first one is that for a ^— dimensional quadratic recursion 
formula, the number of coefficients to be determined at 
order m grows like Z™ and consequently optimization is 
an important consideration. Second, the problem has no 
solution if a given eigenvalue can be written exactly as 
a product of other eigenvalues. This is the exponenti- 
ated form of the problem of logarithmic anomalies raised 
by Wegner in Ref. Generically, such a problem is 
likely to occur is an approximate way. For instance, if 
we use the eigenvalues for the D = 3 hierarchical model 
given in Ref. we find that A3 - A2 ~ 10"^. When 
this is the case, we have a "small denominator problem" 
which refiects approximate "resonance" among the vari- 
ous "modes" present. Pursuing this analogy, the results 
presented here provide a solution of a non-linear problem 
with one degree of freedom. Their application to realistic 
systems seems likely to have a complexity and an interest 
comparable to systems of coupled non-linear oscillators. 
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